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Abstract 

We study the effects of dissipation or leakage on the time evolution of Grover's algorithm for 
a quantum computer. We introduce an effective two-level model with dissipation and random- 
ness (imperfections), which is based upon the idea that ideal Grover's algorithm operates in a 
2-dimensional Hilbert space. The simulation results of this model and Grover's algorithm with 
imperfections are compared, and it is found that they are in good agreement for appropriately 
tuned parameters. It turns out that the main features of Grover's algorithm with imperfections 
can be understood in terms of two basic mechanisms, namely, a diffusion of probability density into 
the full Hilbert space and a stochastic rotation within the original 2-dimensional Hilbert space. 
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Recently, quantum computing has emerged as one of the most challenging fields of physics 
both for theoreticians and experimentalists (see Ref. [1] for a review). At the core of the 
theoretical side, a few quantum algorithms are now available, which can solve a certain class 
of problems faster than any available classical counterparts: for example, Shor's algorithm! 
factorizes a given large number AT at ~ (logA^) 2 time steps with an exponential speed-up. 
Using Grover's algorithm (GA),i one can find a specific item on a long list of size A" at 
~ \/N time steps, which is a considerable gain in speed as compared with ~ A" in classical 
algorithms. 

These quantum algorithms operate perfectly only on ideal quantum computers. On the 
other hand, a certain amount of dissipation or uncontrolled coupling to the environment 
is clearly inevitable on real quantum computers. For example, any deviation from ideal 
operation in quantum gates, which may result from various origins, including fluctuation in 
the excitation energies of two- level systems (qubits), can be considered as "imperfections". 
The imperfections will affect the efficiency of a quantum computer, and the operability of 
a given quantum algorithm may break down to the point of losing its advantage over a 
classical counterpart. Therefore, it is of vital importance to have a sound picture of how 
an error due to the presence of imperfections evolves in quantum algorithms. Obviously, a 
reasonable picture of the basic mechanisms given by the imperfections will be very crucial 
in constructing an appropriate quantum error correction method. 

In general, the quantum state in a quantum computer is essentially a many-body (or 
network) state, the time evolution of which is delicately controlled by a given quantum algo- 
rithm. From such a point of view, the study of imperfection effects on quantum algorithms 
would belong to a more general research field which investigates disorder effects on the dy- 
namics of a many-body state. Their exact treatment is actually a complicated subject, and 
only a few results have been obtained giving either general frameworks for understanding 
the effects or general methodologies for calculation. 

There exist several theoretical, mainly numerical, investigations in this direction. The 
main stress has been given, from a practical point of view, on the stability of quantum 
algorithms with respect to the presence of imperfections. Cirac and ZollerS reported 
that the operability of quantum computing is rather safe against disorders available in 
the quantum Fourier transform process. In Refs. [5] and [TU], the disorder effect in Shor's 
algorithm applied to the factorization of the number 15 was studied and by using the fidelity 
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being denned as the square of the overlap of the actual quantum state with the ideal one, 
it was found that the operability of the Shor's algorithm can be destroyed due to a very 
small strength of the disorder in the modular exponentiation part.0 More systematic results 
have recently been obtained in Ref. |TT| from the study of quantum computing of quantum 
chaos and imperfection effects: by considering the presence of imperfections in the quantum 
Fourier transform, it was obtained that the imperfection strength scales polynomially with 
the number of qubits for the inverse participation ratio (IPR), which measures the strength 



of localization of quantum state and plays a role of the fidelity in Ref. Nevertheless, it 
still remains at a primitive stage regarding an understanding of basic mechanisms carried 
by the imperfections in quantum algorithms. So far, the main policy has been simply to 
watch a deviation of the quantum state from the ideal one and to analyze its parameter 
dependence. 

In this paper, we investigate the time evolution of a state governed by Grover's algo- 
rithm with imperfections, with a main emphasis on an understanding of interplay of the 
imperfections with the algorithm operator. Based on the idea that the ideal GA operates in 
an effective 2-dimensional Hilbert space, a stochastic two-level model with dissipation will 
be introduced, and then its simulation results will be compared to those of the GA with 
imperfections, which operates in a larger relevant Hilbert space resulting from the presence 
of the imperfections. They are in a good agreement via an appropriate fit of parameters. 
An analytic solution of the two-level model is given with some modification and provides a 
comprehensive picture of imperfection effects on the GA. 

Let us begin with a brief sketch of the GA. The final goal is to identify \j) (target state) 
among N = 2 nq quantum states, where n q is the number of qubits. Initially, the state of 
quantum register is prepared as a superposition of all states with the same amplitude. The 
GA may be broken up into two steps: (i) rotation of phase of \j) by 7r and (ii) application 
of a diffusion operator D which is defined, in matrix form, as Dki = —Ski + 2/N with 
k, I — 0, 1, • • • , N — 1, and 5m denoting the Kronecker delta. The step (ii) is achieved by 
applying the Hadamard operation to each single qubit and then performing a conditional 
phase shift on the computer with every computational basis state except \k — 0) receiving 
a phase shift of —1 followed by the second Hadamard operation to each single qubit. Then, 
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the quantum state during time evolution can be expressed a 
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|¥(0)> = sin^b) + ^==J2\k). (1) 



The initial state is characterized by § = $o with sin^ = 1/yN. Each iteration transforms 



into + a;)), where sincu = 2a/^V - 1/N. Then, after m w (7r/4)ViV iterations, -8 
becomes very close to tt/2, and a measurement of the state yields \j) with an error 0(1/N). 
We note that the evolution of according to the GA is restricted to a 2-dimensional 



Hilbert space which is spanned by \x) = (1/y/N — 1) J2k^j \k) and \y) = Each iteration 
represents a rotation of the quantum state by the angle u in the x-y plane and the Grover's 
operator for a single iteration can be written in a familiar form 

* , s I cos u — sin u \ 

R(u) = \ (2) 
I sin uj cos uj J 

on the basis {\x), \y)}. 

Imperfections are introduced in the GA as follows: the ideal Hadamard operator in the 
step (ii) is given by n ■ a, where n = (l/y/2, 0, l/v2), and (? x { y ,z) denotes the Pauli spin 
matrix. We now replace n by 

fh q = —j= (cos ip q ■ sin Sg + cos 8 q , a/2 sin <p q ■ sin 8 q , — cos ip q ■ sin S q + cos 5 q ) , (3) 

where q = l,2,---,n q represent each single qubit. Here, 5 q and ip q with \5 q \ < e/2 and 
< (p q < 2it are randomly chosen in an iteration of the GA and also vary randomly from 
iteration to iteration. Then, it turns out that m q is a unit vector tilted from n by ~ e. 
It should be noted that, in spite of the imperfections, since the quantum state evolves 
without coupling to the additional environment, the qubit rotations remain unitary, keeping 
the normalization condition (^1^) = 1 for any iteration number t. The presence of the 
imperfections will provide an additional coupling between the 2-dimensional Hilbert space 
spanned by {\x), \y)} ("computational space") and the rest part of the total Hilbert space 
with 2 nq dimensions, leading to the quantum leakag^ from the computational space as an 
intrinsic source of error in ideal gate operations. 

Typical results of the GA with imperfections are shown in Fig. 1: (pj) and F denote an 
ensemble- averaged probability of the target state \j) and an ensemble-averaged fidelity over 
100 random runs, respectively, each of which is here given for n q = 13 and for imperfection 



strengths e = 0, 0.005, 0.01 and 0.02, respectively, as a function of iteration number t. 
Clearly, they are given by 



respectively, where the outer bracket represents the ensemble average. In the case of e = 0, 
{Pj)o(t) oscillates between and 1 and reaches 1 at t (m+l/2)(n /2)\/N 71, 213, 355, • • • 
with m — 0, 1, 2, • • •. When e is non-zero, one still finds oscillating features with the same 
period as in the ideal case, however, with an envelope decaying nearly exponentially with 
time t. As t increases, the system approaches a saturated regime, where the noise completely 
dominates the ideal system dynamics, and accordingly (pj) e (t) fluctuates around 1/N. A 
novel feature is that the decay affects not only the shape of the upper envelope but also 
that of the lower envelope such that the lower envelope is not simply given by (pj) € (t) = 0. 
This means that the probability for the system to remain at the target state is still available 
even at the time it originally vanishes in the ideal unperturbed system. Furthermore, F e (t) 
is found to approximately equal the upper envelope of {pj) e (t). 

As noted earlier, in the absence of imperfections, the wave-function of quantum register 
evolves within a very small part (of dimension 2) of the total Hilbert space (of dimension 
2 riq ). Furthermore, since the amplitudes of \x) and \y) remain real or at least keep the 
same phase over the time evolution, the actual relevant space is even smaller than the 2- 
dimensional entire Hilbert space. Let us denote the 2-dimensional Hilbert space spanned 
by \x) and \y) and the total Hilbert space by 7i 2 and 7i t , respectively. The above results 
suggest that in general, the disordered GA operator yields states which are not restricted 
in H2 but spread over a larger space Tit ("computational leakage"). In other words, the 
presence of imperfections induces a probability density flow from Tt 2 to Tit with diffusion- 
like nature. Then, let us define |u> 2 (t)| 2 as the probability that the state remains in Tt 2 at 
time t with an exponentially decaying function of t, 



where 7 represents the strength of the diffusion which depends on system parameters such 
as the strength of imperfections and the qubit numbers. Also, the imperfections affect the 
dynamics of the state within Ti 2 - in general, the phases of the two amplitudes of \x) and 
\y) are not equal to each other, and it is reasonable to assume that random phases are 



(Pj)e(t) = (|<J|*M)>| 2 ) , F £ (t) = (|<*(e = 0,t)|*(e,t)>| 2 ), 



(4) 



w 2 (t)=e^\ 



(5) 



5 



introduced during each iteration. Therefore, we would now like to adopt an effective two- 
level model which can encapsulate the effects of imperfections in the GA living in 7i t - Here, 
the time evolution of a quantum state \il>(t)) = c m (t) \m) + c n (t) \n) on the bais {\m), \n)} 
is described by 



c m {t + l) 



= e 7 R(u) ?7(0 m ,0 n ) 

^ c n (t + i) j y c n (t) 



(6) 



where U((p m , 4> n ) is a diagonal matrix with U mm = e l ^ m and U nn = e^ n , and 4> m (t) and 4> n if) 
are assumed to be two independent random variables without any time correlation. Let each 
of these phase variables be chosen from a box distribution [— W^/2, W^/2] for a given W^. 
The frequency uj = sia~ 1 (2y/N — l/N) is the same as in the GA, and the initial conditions 
are given by c m (0) = cos^o and c n (0) = sin^o with $0 = sin _1 (l/ yN). This is a stochastic 
two-level model with dissipation and we refer to it as STLM hereafter. Here, we obtain, 
after a minor calculation, an ensemble- averaged probability of the target state \j) and an 
ensemble- averaged fidelity, respectively: 

(Pj)w\(t) = (\c n (t)\ 2 ) , F™(t) = (\c m (t) cosM + tfo) + c n (t) sin(a;t + ^ )| 2 ) . (7) 

It is noteworthy to consider the difference between the STLM and the original GA with 
imperfections. First, the finite fraction 2 _n « +1 occupied by Ti.2 in 7i t is neglected in the 



STLM so that w%(t) decays to zero instead of ~ y / 2~ n i +1 . Since we are interested mainly 
in the regime before saturation, this is clearly not a significant difference. Secondly, the 
stochastic features of 7 are not considered. But, this is not critical, either, since those 
features will contribute a negligible correction to 7 after an ensemble-average in eq. (0). 

Now, we perform a numerical simulation to obtain (pj)w\(t) an d Fyfi(t), which will be 
compared with (pj) e (t) and F e (t) of the GA with imperfections, respectively. In Fig. 1, 
the results from the STLM are shown as solid lines: they are given by ensemble-averages 
over 1000 realizations, respectively. We find that these results from the STLM provide an 
impressive agreement with the results of the GA after a proper adjustment of 7 and W^. This 
suggests that the main physical ingredients of the disordered GA are correctly incorporated 
in the STLM. Nevertheless, the origin of the novel feature in the lower envelopes is still 
unclear. 

Without loss of generality, (c m (t), c n (t)) in the STLM during the time evolution can be 
written by (e~ 7 * cos^t), e~" /t+l ^ sin^t)) with (p(t) := <p n — <p m up to an overall phase. In 



case of <j)(t) = for arbitrary t, the angle $(t) increases by u after each iteration and is then 
given just by ut + tV However, if </>(t) does not vanish, then from its stochastic nature, it 
follows that — — 1) is not constant but would fluctuate around uM Now, under the 
assumption that and <p(t) are not correlated with each other, but simply two random 
variables, we can find analytic expressions of (Pj)w(t) and F${t), respectively; let <&{t) 
increase by u + rjt-i between t — 1 and £ such that 

t-i 

= tf + + E Vk, (8) 

fc=0 

where % forms a GAussian distribution with mean and width A$, and then Y^k=oVk also 
satisfies a GAussian distribution with mean and width A^y/t. From this and eq. (0), we 
get: 

(p.-)W(t) = |u> 2 (t)| 2 (sin 2 ^(t)\ = |W2 ^ /°° sin 2 (^ + ^ + x) e"" 2 /^') rfx 

X 7 A^y/nt J-oo 

= — — [1 -cos(2cjt + 2^ ) -e- A ^\ (9) 

(note that no subindex W<p appears in (pj)^(t) ). If we further assume that <f>(t) also is of 
a GAussian distribution with mean and width A^, we then arrive at 

F™(t) = "—^ [l + e- A ^{l-sin 2 (2^t + 2^ )-(l-e- A ' /4 )}] . (10) 

Fig. 2 shows a comparison between the results of the GA with imperfections and those of 
eqs. (|9p and (0). The good agreement in (pj)(t) would provide an explanation of why its 
lower envelope is not simply given by (pj) = in the GA with imperfections; the uncertainty 
in the rotation angle during a single iteration accumulates as the iteration proceeds. Then, 
does not represent a definite direction on a 2-dimensional plane but spreads over an 
interval range (—A$y/i, A^y/i). This offers an additional decay channel into the target state 
\j) after ensemble-averaging (see the term e _A2,< in eq. (f|)). Also, in eq. flT(i|) with = 
we have FW(t) = (e" 2 ^/2) (l + e A 5» , which is given by the solid lines in Fig. 2 as the 
best fit of F e (t) of the GA with imperfections. From eq. (|I0|) , it immediately follows that 
F(A ^ 0) is always less than F(A^ = 0). 

In summary, we have investigated imperfection effects on the time evolution of the 
Grover's algorithm both numerically and analytically. An effective two-level model with 
dissipation and randomness has been introduced and the results show a good agreement with 



the simulation results of the disordered Grover's algorithm. It turns out that the main fea- 
tures in the results of the disordered Grover's algorithm can be understood through the 
diffusion-like behavior of quantum states from the original partial Hilbert space into the 
full Hilbert space. The two main decaying mechanisms found in this work are its direct 
manifestations. Our finding will provide a useful basis for study of more general imperfection 
effects in quantum algorithms. 
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FIG. 1: Behaviors of (pj) (•) and F (o) in the Grover's algorithm for the qubit number n q = 13 
with the imperfection strength (a) e = 0.005, (b) e = 0.01 and e = 0.02 (inset). Each data point 
is given, for every 20 iterations, by an ensemble-average over 100 realizations. The dotted line of 
(a) represents (pj) for n q = 13 in the ideal case (e = 0). The solid lines result from the stochastic 
two-level model described in the text with parameters (a) 7 = 7.6 x 10 -4 , (b) 3.0 x 10 -3 and 
1.3 x 10" 2 (inset), and (a) = 0.089, (b) 0.19 and 0.25 (inset), respectively. 
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FIG. 2: Comparison between the results of the Graver's algorithm with imperfections and the 
theoretical predictions given by eqs. @ and ([!(]). The symbols indicate the same data as in Fig. 1. 
The parameters 7's are the same as in Fig. 1 with (a) A$ = 2.0 x 10~ 2 , (b) 4.2 x 1CT 2 and 3.5 x 1CT 2 
(inset), respectively, and = for all three. 
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